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Since  many  methods  of  quantitative  chemical  analysis  are  not  "fully  selective" 
(1) ,  or,  in  other  words,  not  free  from  interferences,  important  goals  of  the 
analyst  are  to  detect,  characterize  and,  hopefully,  eliminate  interferences  that 
arise  during  the  course  of  an  analysis.  Interferences  that  do  not  arise  from 
chemical  species  (e.g.,  electronic  drift)  in  the  sample  clearly  should  be  eliminated. 
However,  interferences  (i.e.,  chemical,  physical  or  spectral)  arising  from  chemical 
species  in  the  sample  can  be  regarded  as  "misplaced  analytical  signal"  and  it  would 
be  best  to  sort  these  interferences  during  a  multianalyte  analysis.  For  example, 
it  is  well  known  that  when  the  422.7  nm  line  is  used  for  the  atomic  emission  analysis 
of  calcium,  the  presence  of  sodium  in  the  sample  will  cause  an  overestimation  for 
the  quantity  of  analyte  in  the  sample  due  to  an  interference.  It  is  common  practice 
to  apply  a  correction  to  the  calcium  signal  based  on  the  strength  of  the  sodium 
signal.  However,  this  negative  correction  to  the  calcium  signal  ideally  should  be 
accompanied  by  a  positive  correction  to  the  sodium  signal.  In  this  way,  the  in¬ 
fluence  of  sodium  on  the  calcium  line  would  not  be  wasted  but  would  be  correctly 
sorted  giving  rise  to  a  conservation  of  analytical  signal  and  increased  sensitivity. 

Ideally,  a  method  is  clearly  needed  that  can  1)  detect  interferences  during 
a  multianalyte  analysis,  2)  mathematically  characterize  these  interferences,  no 
matter  how  complex,  and  3)  sort  all  interference  signals  to  al~ow  the  maximum 
utilization  of  all  analytical  signals.  This  method  would  not  only  provide 
more  accurate  multianalyte  analyses  hut  would  also  serve  to  characterize  the 

t 

selectivity  of  the  analytical  Instrument  and  therefore  be  an  important  aid  for 
analytical  method  development. 

Recently,  our  laboratory  developed  the  mathematical  basis  for  just  such  a 
method.  It  is  based  on  the  preferred  method  of  standard  additions  and  is  a  gen¬ 
eralization  to  Include  any  number  of  analytes.  Hence  it  is  referred  to  as  the 
Generalized  Standard  Addition  Method  (GSAM)  (2)  for  simplicity.  The  GSAM  is 
applicable  to  multicomponent  analysis  where  Interferences  of  any  kind  are  a  problem. 
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This  paper  deals  with  the  limitations  of  the  linear  model  GSAM  and  the  devel¬ 
opment  and  testing  of  optimal  computer  algorithms  and  experimental  designs  for 
carrying  out  multicomponent  analyses  using  the  GSAM.  During  the  course  of  the 
study,  certain  fundamental  concepts  and  limitations  to  multicomponent  analysis 
were  discovered.  These  concepts  allow  the  analyst  to  calculate,  in  advance,  the 
penalty  or  cost  of  performing  a  multicomponent  analysis  when  interanalyte  inter¬ 
ferences  are  present.  This  penalty  takes  the  form  of  an  amplification  of  the 
measurement  error  as  it  is  propagated  to  the  final  results  (initial  analyte  con¬ 
centrations  and/or  response  constants).  This  theory  should  be  of  considerable 
utility  to  the  analytical  chemist  as  it  describes,  quantitatively,  the  usefulness 
of  an  analytical  system  and  can  even  be  used  to  compare  and  select  systems.  This 
same  theory  has  vastly  improved  our  understanding  of  the  characteristics  of  the 
GSAM  and  has  provided  direction  for  the  optimization  of  its  application. 

In  the  remainder  of  this  paper,  the  linear  response  GSAM  is  briefly  reviewed, 
the  theory  mentioned  above  is  described  in  detail,  and  the  improvements  given  by 
the  theory  are  tested  on  a  four  component  spectrophotometric  analysis. 

BACKGROUND  AND  THEORY 

The  Generalized  Standard  Addition  Method  -  A  Brief  Overview 

As  any  analyst  knows,  the  standard  addition  method  requires  the  measurement  of 
the  analytical  signal  before  and  after  addition  of  known  quantities  of  the  analyte 
(standard  additions)  have  been  added  to  a  sample.  The  analyte  concentration  is 
found  by  extrapolating  the  regression  line  to  the  negative  portion  of  the  plot  of 
signal  vs.  concentration  added  (Ac).  The  model  assumes  that  the  analytical  signal, 
or  more  generally  the  response,  R^  is  given  by: 

R^  -  ck^  -  (Ac  +  c, )  k^  1) 

where  Ac  is  the  effective  total  concentration  change  of  the  analyte  after  any 
standard  addition,  c«  is  the  initial  analyte  concentration  and  k^  is  the  linear 
response  constant.  Of  course,  if  analytical  signal,  H,  is  responsive,  again  in  a 


linear  manner,  to  r  components  in  the  sample  (interference  effects),  the  true 
response  equation  is: 


h  -  l< =sk.l  J> 

s-1 

where  c  is  the  concentration  of  the  s— -  component  and  Its  response 
constant.  For  later  use,  it  is  most  important  to  note  that  can  be  expanded 
as: 

h  ~  *  2.  •csk,*  3) 

S*1  S*1 

which  simply  says  that  the  total  concentration  is  the  sum  of  each  initial  concen¬ 
tration  plus  its  change  in  concentration  brought  about  by  standard  additions. 

Now,  for  r  analytes,  p  analytical  responses  (p  r)  must  be  measured.  Each 
response  can  depend  on  any  combination  of  the  r  components  provided  that  each  of 
the  r  components  (hereafter  referred  to  as  analytes  as  all  r  concentrations  will 
be  determined  simultaneously)  affects  at  least  one  of  the  p  responses.  The  GSAM 
requires  that  n  multiple  standard  additions  (MSA's)  be  made  so  as  to  span  the 
r-dimensional  concentration  space  (n  _>  r) .  Expressing  the  above  in  matrix  form 
gives , 

4) 

where.fi,  is  the  n  x  p  matrix  of  measured  responses  to  the  n  MSA’s.  is  the  n  x  r 
matrix  with  the  total  concentration  (Ac  +  c0)  of  the  r  analytes  in  each  of  the  n 
rows  andj£  is  the  r  x  p  matrix  of  constants  relating  the  contribution  of  each  of 
the  r  analytes  to  the  p  responses.  Note  that  many  elements  of  £  can  be  near  zero 
If  interferences  are  not  present. 

The  solution  for  is  given  in 
reference  (2)  to  be, 

(ACT  AC)-1  ^£TfiR,  5) 

where  ££  is  an  n  x  r  matrix  of  total  concentration  changes  for  r  analytes 


at  each  of  n  MSA's.  Given  c..  the  column  vector  of  initial  analyte  concentra¬ 
tions,  can  be  found  by 

c/  -  r^1  jC1,  r*  -  (.Rj,  •••,  .Rr>T  6) 

It  is  very  important  to  recognize  that  AC  in  equation  5  is  the  matrix  of 
effective  total  concentration  changes  for  each  analyte  after  the  MSA’s  and,  unless 
volume  changes  are  negligible,  is  cannot  be  known  since  the  initial  concentrations 
are  unknown.  In  practice,  this  is  not  a  problem  as  simple  volume  corrections 
convert  concentrations  to  absolute  quantities  (e.g.,  moles). 

Equation  2  can  be  rewritten, 

K,  ■  l  ,  7> 

*  s-1  V 

where  N  is  the  total  number  of  moles  of  analyte  s  in  a  volume  V.  This  leads  to 


«» 5  *  I,  8> 

s=*l 

which  can  now  be  separated. 

On  -  +  l  ’\\l  ’> 

8=1  S=1 

0N  is  the  initial  number  of  moles  for  analyte  s  and  AN  in  the  total  amount  of 
s  s 

analyte  s  added  which  is  now  known.  Hence,  for  equation  4,  if  all  responses  at  all 
MSA's  are  volume  corrected, £  is  the  matrix  of  amounts  expressed  in  moles  and  the 

solution  can  proceed  as  in  reference  2.  The  solution of  the  overdetermined  system 

10) 

is  now 

T  -1  T  11) 

K  -  (AN  AN)  1  AN  A^ 

with  AN  expressed  in  total  changes  of  quantitites  and  A^ equal  to  the  matrix  of 
volume  corrected  response  changes. 

The  GSAM  requires  that  MSA's  be  made  prior  to  each  of  the  n  measurement  steps. 
One  MSA  may  consist  of  the  addition  of  one  standard  or  several  standards.  The  MSA's 
are  attempts  to  span  an  r-dimensional  concentration 


space.  If  the  linear  model  holds  for  the  simple  standard  addition  method  (one 
analyte)  the  method  of  additions  is  clear  (3).  Unfortunately,  this  is  not  the 
case  for  the  GSAM. 
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The  goal  in  this  study  was  to  determine  the  optimal  method  of  performing  the 
multiple  standard  additions.  "Optimal"  meaning  the  lowest  errors  in  the  determina¬ 
tion  of  K  and  n^. 

The  accuracy  of  the  determination  of  no,  the  vector  of  initial  amounts  of  all 
the  analytes,  is  dependent  on  at  least  five  items. 

1)  The  accuracy  of  each  response  is  of  paramount  importance  and  represents  a 
limiting  condition  as  in  any  quantitative  analysis. 

2)  With  the  GSAM,  the  multiple  standard  additions  must  be  made  without  deter¬ 
minate  error  and  with  acceptable  precision.  Again,  this  is  true  of  any  analysis 
regardless  if  calibration  or  standard  addition  methods  are  used. 

3)  The  third  item  deals  with  interferences  which  can  be  described  using  math¬ 
ematical  properties  of  the  K  matrix  to  be  explained  later. 

4)  The  experimental  design,  as  might  be  expected,  can  strongly  influence  the 
accuracy  of  K  and  hence  n^ 

5)  The  mathematical  algorithms  selected  for  computation  can  effect  the  results 

due  to  round-off  errors.  This  will  be  clear  in  the  following. 

The  emphasis  of  this  paper  is  item  4.  However,  in  a  search  for  better  ways  to 

implement  the  GSAM,  much  has  been  learned  about  all  five  items  and  a  more  fundamental 
understanding  of  multicomponent  analysis  has  been  the  result.  This  paper  describes 
the  manner  in  which  measurement  accuracy  (item  1),  experimental  manipulation  involved 

with  calibration  (item  2)  and  interferences  (item  3)  Influence  chemical  analyses, 
and  how  the  way  in  which  the  analyses  are  performed  (item  4)  can  alter  these 
Influences  and  the  results  calculated  (item  5) . 

Error  Propagation 

After  computation  of  the  K  matrix,  the  desired  initial  quantities  (the  com¬ 
ponents  of  no)  are  obtained  by  solving 

T 


* 


«i«Kll'W>  tn  i 


- - - _____ - 
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If  the  matrix  Is  well  conditioned  (i.e.,  the  analytical  sensors  are  very 

selective)  any  linear  equation  solver  may  be  applied  to  equation  12.  If  not 
(i.e.,  the  columns  of  £  are  nearly  linearly  dependent),  Househoulder ' s  triangular 
decomposition  (see  Appendix)  should  be  preferred  in  order  to  at  least  guarantee  a 
minimum  of  numerical  rounding  errors  (item  5  above).  Nevertheless,  the  magnifica¬ 
tion  of  experimental  errors  may  be  crucial. 

The  experimental  application  of  the  GSAM  is  as  subject  to  experimental  errors 
as  any  measurement  procedure.  These  errors  may  be  propagated  twice:  first, 
in  determining  the  K  matrix  by  solving  the  least  squares  problem  and,  second,  in 
solving  the  linear  system  12  using  the  calculated  K  matrix.  The  sensitivity  of  the 
solution  of  a  numerical  problem  with  respect  to  errors  (i.e.,  the  size  of  the 
error  magnification  factor)  is  usually  called  the  condition  of  the  problem.  Even 
in  the  linear  case  the  relative  data  error  may  be  magnified  by  a  factor  which 
can  lead  to  completely  unreasonable  results.  Hence,  a  very  careful  analysis  should 
be  applied  to  the  sensitivity  of  the  final  result  n0  with  respect  to  instrumental 
and  performance  errors  and,  consequently,  to  the  optimal  design  of  the  addition 
matrix  AN. 


Error  amplification  can  be  shown  by  a  simple  but  illustrative  example.  Assume 
the  following  two  by  two  system  of  linear  equations,  Ax  *  b,  is  to  be  solved. 


T 

The  solution  is  x  *  (-4.0,  3.0)  as  is  easily  verified.  Now  suppose  the  right-hand 

T 

side,  ^  is  perturbed  by  the  error  vector  (0.3,  0.4)  which  produces  a  relative 

error  ||^||  /||£||  of  2%.  (||  •  ||  =  square  root  of  sum  of  squares  of  vector  elements.) 

The  corresponding  solution^  +  t0  ^  +  *s  A  “  (-3.7,  3.4)T  which 

corresponds  to  a  relative  solution  error  ||jy^||  /||jJ|  of  10%.  Thus  an  error  amplifi¬ 
cation  factor  of  5  results.  To  put  this  result  in  a  more  general  context  the  intro¬ 
duction  of  the  norm  of  a  matrix  is  required. 


mm 


For  a  matrix  AeK  *  ,  the  number 

1 1  ll2 

1 1  i  1 12  -  in*3  ‘  “  ®ax  ll^illo  (lub  means  l_owest  upper  bound) 

iemr  IUII2  HiA"1  '  ‘ 

will  be  called  the  Euclidian  norm  (or  simply  norm)  of  A.  The  square  roots  of 


the  eigenvalues 


T 

of  the  positive  definite  matrix  i  £e!Rr  ,r  are 


called  singular  values  of  In  the  case  where  n-r  and  A  is  normal  (i.e., 

T  T 

^  A  *  ^A  )  the  singular  values  are  the  absolute  eigenvalues  of  i  (e.g.,  all 
symmetric  matrices  are  normal) .  It  is  well  known  that  the  biggest  and  the  re¬ 
ciprocal  smallest  singular  values  of  £  are  equal  to  the  norms  of  A  and  A-1, 
respectively  (4). 

II  A||-  ll  A_1  H  “  °n~1-  1- 

Using  the  same  notation  as  before  and  starting  with  equation  12  where  K,  the 
matrix  of  response  constants,  is  assumed  known  exactly  for  now,  and  q0  is  the 
vector  of  initial  volume  corrected  responses,  a  useful  result  can  be  found. 

It  is  known  that 

llj.ll  £  II i’ll  •  Ha.il  it 

If  1  exists,  as  it  must  in  chemical  analysis,  then 
T  -1 

n_&  “  (K  )  q,  1 

Now,  a  small  change  in  q0,  represented  by  6qe,  gives  rise  to  an  error  6n-, 
which  must  be  kept  small  in  a  chemical  analysis.  Since, 

j5ng  -  QCT)  1  6q,,  1 

II  6n.  II  i  II  (kV1!!  •  II  8a-l(  1 

Combining  16  and  19  and  using  ||  Q^)~*||  -  ||£*||  give 

11  „  „  „  _i„  liia-li 

— — —  1 II X  II  •  Hi  II  •-==-  2 

II  n.  II  II  q.  II 


on,  11  , 

T~-  £  Hi  II  •  Hi1! 

II  n*  II 


Defining  the  condition  number  of  any  nonsingular  matrix  A  as 
cond  <£)  -  ||i||  •  ||i_1 1| 


21) 


8 


and  using  Che  fact  Chat 


ll  5n.ll 

-  <  cond  (K) 

Hn.lL" 


T 

cond  (A)  -  cond  (A  )  the  final  result  Is 

II  gqoll 
II  q.  II 


22) 


The  condition  number  can  have  different  values  depending  upon  the  chosen  norm.  It  is 
always  greater  than  one  and  herein  lies  an  important  result.  If  q.  is  the  true 
response  vector  that  gives  rise  to  the  true  analyte  vector  n0,  then  22 
shows  that  a  small  error  in  response  6q.  can  be  magnified  by  cond  (K^  in  the  worst 
case  to  produce  a  greater  relative  concentration  error;  error  amplification. 

All  of  the  above  relates  to  the  third  item  mentioned  in  the  introduction.  The 
analyst  selects  a  particular  analytical  method  characterized  by  the  matrix  K,  the 
response  constants  of  the  method.  Inequality  22  characterizes  the  "cost  of  inter¬ 
ferences"  in  the  method.  The  cost  is  paid  in  increased  relative  uncertainty  in  the 
initial  concentrations  to  be  determined.  The  ideal  case  of  fully  selective  sensors 
(no  interferences)  with  equal  sensitivity  gives  rise  to  a  £  matrix  with  a  condition 
number  of  one  and  absolutely  no  error  amplification. 

During  an  actual  application  of  the  GSAM,  the  K  matrix  is  determined  as  the 

solution  of  an  overdetermined  linear  system  and  therefore  is  not  known  exactly. 
Hence,  the  propagation  of  errors  due  to  the  K  matrix  must  be  examined.  The 

following  similar  error  estimate  taking  matrix  entry  errors  into  account  can  be 

found  in  reference  4»  p.  195: 


Hi -ill 
Hill 


cond 


1-cond 


23) 


where  x  is  the  solution  of  the  perturbed  system  <^A+£)i  *  whose  error  matrix  £ 
has  to  be  sufficiently  small  (lli^ll  K  D  • 

Using  23  and  22  (by  a  Taylor  expansion  neglecting  higher  order  terms)  the 
combined  error  estimate  is  obtained 


II  5n- 1| 


<  cond 


n, 


© 


II  5a. II  II  5K||" 

Hi*  II  III  II  ^ 


24) 
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Now,  to  estimate  the  second  term  in  the  brackets  a  sensitivity  analysis  of  the 
generalized  solution  of  equation  11  with  respect  to  experimental  error  amplification 
must  be  done.  It  is  considerably  more  complicated,  whence  only  the  final  inequalities 
are  cited  here.  For  details  the  reader  is  referred  to  (4,  5,  8).  The  first  estimate 
analogous  to  equation  22  (transcribed  to  this  paper's  notation)  can  be  found  in 
Reference  4,  p.  221,  to  be, 


k-^ll  llAq.-Aqjl 

-  <  cond  (AN)  — 

ii  1 1  —  'mam' 


where  cond  (^AN)  *  ||  AN||  ||  AN  ||  is  now  the  generalized  condition  number  of  A^  and 

Aq j  and  Aq^  are  the  projections  of  Aq  and  onto  the  range  of  AN.  The  generalized 

condition  of  the  rectangular  matrix  AN  may  be  computed  by  the  formula 

t  1  /2 

cond  (AN)  *  [cond  (AN  AN)]  ,  26) 


where  the  number  in  brackets  denotes  the  usual  (see  equation  21)  condition  of  the 
T 

square  matrix  AN  AN. 

Further  generalizations  are  not  only  beyond  the  scope  of  this  paper  but  are 
actually  unnecessary  for  achieving  our  goals.  The  analytical  chemist  is  not  usually 
interested  in  calculating  upper  bounds  on  the  relative  errors  of  initial  concentra¬ 


tion  vectors.  Actually,  statistical  error  propagation,  the  subject  of  other  work 
in  our  laboratory,  is  usually  of  greater  interest.  Rather,  the  forms  of  equations  24 
and  25  are  far  more  important  as  they  show  the  dependence  of  analytical  accuracy, 

£no,  on  measurement  accuracy,  6qe,  interferences,  cond  and  the  experimental 

design,  cond(^AN).  Functional  dependence  of  item  two  is  covered  in  the  Appendix,  its 
effect  being  a  bit  more  complex.  These  dependences  demand  that  measurement  errors 
be  small,  interferences  be  few  and  small  [cond  close  to  one],  multiple  standard 
additions  be  made  as  accurately  as  possible  (Appendix)  and  that  the  experimental 
design  lead  to  the  smallest  cond  as  possible.  The  first  thr.ee  demands  are  not 

unknown  to  the  analytical  chemist.  But  the  final  demand  has  guided  this  work  to  the 
experimental  designs  proposed  In  the  following  section. 
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Suggested  Variants  of  the  GSAM 

The  sensitivity  of  the  initial  quantities,  n0,  is  highly  dependent  on  the 
condition  of  and,  thus,  on  the  experimental  design  of  the  GSAM.  In  order  to 
minimize  the  error  propagation  factors  two  variants  of  the  GSAM,  as  presented  in 
the  original  paper,  are  now  suggested  for  use.  ■ 

The  original  GSAM  method  used  the  total  difference  calculation  and  does  not 
differ  in  experimental  design  from  the  now  suggested  incremental  difference  calcula¬ 
tion.  Rather,  the  difference  lies  only  in  the  definition  of  and,  correspondingly, 
of  AQ.  The  entries  in  the  AN  matrix  for  the  total  difference  calculation  are  the 
accummulative  amounts  of  standards  added.  For  the  incremental  difference  calcula¬ 
tion,  the  entries  are  just  the  amounts  added  at  any  one  MSA.  Likewise,  the  AQ 
entries  for  the  total  difference  computation  are  the  total  changes  in  response  from 
q,  and  only  the  stepwise  incremental  response  changes  for  the  incremental  difference 
computation. 

The  second  variant  procedure,  called  the  Partition  GSAM,  or  PGSAM,  amounts  to 
a  true  change  in  experimental  design.  The  sample  is  first  split  into  as  many 
partitions  as  there  are  analytes.  (In  practice,  the  analyst  would  begin  with  r 
aliquots  of  the  sample.)  Then,  for  each  partition,  all  of  the  responses  are  measured 


before  and  after  standard  additions  of  only  one  standard  per  partition. 

As  will  be  seen  in  a  later  section,  the  condition  of  ML  for  the  incremental 

difference  computation  GSAM  and  the  PGSAM  is  unity  amounting  to  no  error  amplification 
for  either  variant.  (Actually,  this  is  true  only  when  identical  quantities  of  each 

standard  are  added  or  is  scaled  in  a  manner  to  be  described  later.)  The  normal 
GSAM  is  performed  in  one  vessel  and  the  PGSAM  requires  multiple  vessels.  In  most 
cases,  the  GSAM  would  be  preferred  but  if,  for  example,  additions  of  several  different 
standards  to  a  single  sample  leads  to  changes  in  the  k^'s  (equation  3),  this  model 
nonlinearity  may  be  avoided  by  employing  the  PGSAM.  Experimental  considerations 
can  now  dictate  which  procedure  to  use  and  the  sensitivity  of  the  sensors  can 
dictate  the  amount  of  each  standard  to  use  for  any  one  MSA. 
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EXPERIMENTAL 
Reagents : 

Four  aqueous  solutions  of  0.01  M  NiClj,  0.01  M  CuCl2>  0.01  M  CoCl2,  and 
0.002  M  K2Cr202  (all  A.C.S.  Reagent  Grade)  were  used  to  serve  as  standards  for 
the  additions  and  to  prepare  the  "unknown"  mixture  consisting  of  1.0  ml  each  of 
the  N1C12,  CuC12  and  CoCl2  standard  solutions.  Standard  additions  were  made  in 
units  of  one  ml  of  the  respective  solutions  or  1.0  x  10~^  mole  of  NiCl2>  CuCl2,  CoCl2 
and  2.0  x  10  ®  mole  I^CrjO^. 

Apparatus : 

A  Cary  219  spectrophotometer  was  used  for  all  measurements.  The  responses 
measured  were  the  absorbances  of  every  solution  at  five  wavelengths.  The  wave¬ 
lengths  corresponded  to  absorbance  peak  maxima  for  each  of  the  four  analytes  with 
a  second  peak  maximum  for  nickel  included.  Wavelength  were  assigned  as  follows: 

-  351.0  nm  (Cr^2-) ,  X2  -  394.5  rnn  (Ni2+),  X3  -  511.5  (Co2+),  X4  =  660.0  nm 
(Ni2+),  X5  *  820.0  nm  (Cu2+) . 

Procedure : 

_A:  (GSAM)  16  additions,  4  for  each  analyte,  were  added  to  the  unknown  mixture  in 

the  following  order:  4  times  1  unit  of  NiCl2;  4  times  2  units  CuCl2; 

4  times  4  units  of  CoCl2;  and  4  times  8  units  of  K2Cr202.  After  each 

addition  the  absorbances  at  all  five  wavelengths  were  recorded. 

jJ:  (PGSAM)  Four  equal  aliquots  of  unknown  mixture. 

To  the  first  aliquot  4  additions  of  1  unit  NiCl2> 

to  the  second  aliquot  4  additions  of  1  unit  CuCl2, 

to  the  third  aliquot  4  additions  of  1  unit  CoCl2,  and 

to  the  fourth  aliquot  4  additions  of  1  unit  KjCrjO^,  were  made. 

After  each  addition  the  absorbances  of  the  respective  partition  at  all  five 

wavelengths  were  recorded. 

Computer  Program: 

All  programs  are  written  in  Fortran  IV  and  are  available  from  Infometrix, 

Inc.,  P.  0.  Box  25888,  Seattle,  WA  98125. 
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RESULTS  AND  DISCUSSION 

The  experimental  goal  was  to  corroborate  the  theoretical  considerations 
according  to  possible  addition  techniques.  Hence,  an  experiment  was  selected 
that  was  easy  to  perform  and  that  allowed  the  testing  of  different  addition 
matrices  on  an  instrument  with  a  large  and  linear  dynamic  range.  One  experiment 
that  meets  these  requirements  is  the  determination  of  the  concentrations  of 
inorganic  ions  in  a  mixture  by  visible  light  spectrophotommetry.  As  described  in 
the  Experimental  Section,  two  procedures  were  used  to  perform  the  experiment: 
the  standard  GSAM  and  the  sample  partition  GSAM  or  PGSAM. 

The  results  were  calculated  in  two  different  ways;  the  ^otal  difference 
computation  and  the  iticremental  difference  computation  (TDC  and  IDC) .  During 
the  experiment  the  responses  (absorbances)  at  5^  different  wavelengths  were 
measured  although  only  four  would  have  been  necessary  since  an  unknown  mixture 
consisting  of  only  four  analytes  was  considered.  This  was  done  to  obtain  different 

JK  matrices  using  different  combinations  of  selected  wavelenths  in  order  to 
compare  results.  Of  course,  a  generalized  inverse  solution  for  the  initial  quant¬ 
ities  starting  with  equation  12  can  be  used  when  r  >  p. 

In  a  real  analytical  sample  the  analyst  Is  faced  with  possible  large  differences 
among  the  initial  analyte  concentrations  and  varying  sensitivities  among  the 
transducers  used  to  obtain  the  analytical  responses.  Therefore,  the  concentration 
changes  made  for  each  analyte  may  be  quite  different  in  magnitude.  This  necessary 
element  of  the  experimental  design  produces  large  differences  in  relative  magnitudes 
among  columns  of leading  to  artificially  increased  condition  numbers.  To  get 
the  smallest  possible  condition  number  for  a  certain  addition  matrix  AN  and  thus 
the  smallest  error  estimate  of  thej^matrix  the  addition  and  response  matrices  were 
scaled  such  that, 


l  (An  )2  -  1 

j-1 


fee*  |(l  • 


i  -  1 


27) 


13 


Table  I  shows  the  addition  matrix  and  the  five  volume  corrected  response  columns 

of  the  TDC  before  scaling.  Mote  that  the  additions  are  summed  since  the  total 

change  in  the  response  was  observed  with  respect  to  the  initial  response  Since 

the  additions  are  given  in  above  defined  units,  the  results  are  also  given  in  the 

same  units.  Also,  unequal  amounts  of  each  standard  were  used  for  the  additions. 

Although  this  experimental  design  may  not  be  ideal  from  the  experimental  or 

statistical  point  of  view,  it  gives  a  smaller  cond  (AN)  for  the  total  difference 

GSAM  than  equal  additions,  a  fact  that  easily  can  be  verified  by  the  reader. 

Table  II  gives  the  data  for  the  incremental  difference  computation.  Note  that 

the  addition  matrix  now  consists  of  only  the  additions  made  at  one  time  and  the 

response  matrix  consists  of  only  the  changes  of  the  responses  with  respect  to  the 

T 

response  before  each  addition  was  made.  The  scaled  ^N  matrix  is  therefore  a 

diagonal  matrix  with  condition  number  one. 

Likewise,  Table  111  shows  the  results  for  the  sample  partition  method,  PGSAM. 

The  addition  matrix,  jAl^  shows  that  amounts  are  again  accumulative  but  only  within 

a  single  sample  partition  per  column.  Besides  the  experimental  advantage  mentioned 

T 

above,  the  use  of  the  PGSAM  has  the  additional  advantage  that  the  scaled  AN  AN 
has  a  condition  number  of  one.  The  major  disadvantages  of  the  PGSAM  are  the 
extra  amount  of  sample  required  and  the  manipulation  necessary  when  a  sample  is 
partitioned  which  can  introduce  additional  experimental  error.  In  the  present 
experiments,  the  normal  GSAM  results  were  obtained  using  only  one  pair  of  cuvettes 
for  the  spectrophotometer.  One  reference  cuvette  and  four  sample  cuvettes,  one 
for  each  partition,  were  used  for  the  PGSAM.  Thus  the  errors  associated  with 
volume  transfers  and  cuvette  positioning  in  the  instrument  were  compounded  by  the 
use  of  more  experimental  equipment. 


Formally,  the  first  calculation  step  in  the  GSAM  is  the  determination  of  jC. 
Since  results  for  the  TDC-GSAM,  IDC-GSAM,  and  the  PGSAM  were  nearly  Identical,  onlj 
the  J?  matrix  for  the  IDC-GSAM  for  each  of  the  two  selections  of  wavelengths  is 
given.  Table  IV  shows  the  experimentally  determined  for  the  "most  selective" 
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combination  of  wavelengths.  Although  none  of  the  wavelengths  are  fully  selective 
for  any  one  element,  only  one  wavelength  (Xj)  represents  a  serious  interference 
problem.  The  resulting  error  amplification  limit  [cond  (£)]  is  slightly  less  than 
a  factor  of  two.  Table  V  shows  the  JC  for  a  less  selective  combination  of  wavelengths 
reflecting  more  sericus  interferences  and  a  much  larger  error  amplification  factor 
[cond  (K) ] . 

Table  VI  gives  the  final  results  for  the  TDC-GSAM,  IDC-GSAM,  and  the  PGSAM 
for  the  most  selective  combination  of  wavelengths  and  only  the  IDC-GSAM  for  the 
less  selective  combination  of  wavelengths.  The  result  for  the  second  group  of 
wavelengths  paralleled  the  results  of  the  first  group  and  are  not  reported  here 
for  brevity.  The  eigenvalues  and  condition  numbers  are  all  as  expected.  The  first 
three  estimates  of  cond  00  are  essentially  equivalent  reflecting  only  experimental 
error  of  the  estimates  of  the  K  matrices  pointed  out  above.  The  last  cond  is 
larger  and  quantitatively  estimates  the  potential  cost  of  performing  the  chemical 
analysis  with  a  less  selective  analytical  system. 

The  interesting  results  for  the  analytical  chemist  come  at  the  bottom  of 
Table  VI.  The  most  accurate  results  using  the  first  group  of  wavelengths  is 
achieved  with  the  application  of  the  IDC-GSAM.  Its  average  absolute  %  error  of 
1.33%  is  quite  close  to  the  precision  limit  (est.  1%)  of  the  instrument  used  for  the 
range  of  absorbances  measured.  This  result  also  points  out  a  strength  of  the  GSAM 
that  can  be  understood  by  a  closer  examination  of  Table  II.  Since  the^g  values 
are  volume  corrected,  a  group  of  four  values  in  a  column  corresponding  to  the 
addition  of  equal  amounts  of  any  one  standard  should,  ideally,  be  equal.  Thev  are 
not  equal  because  measurement  and  addition  errors,  in  some  cases,  vary  widely. 

Using  as  true  values  the  averages  of  the  four  sets  of  values  corresponding  to  the 
changes  in  the  primary  spectral  peak  for  each  of  the  four  analytes,  statistical 
standard  error  estimates  of  4.3%,  2.6%,  7.1%  and  3.5%  are  calculated.  These 
standard  error  estimates  are  much  greater  than  the  standard  error  estimates  in  the 
final  concentration  values  (bottom  line  of  Table  VI)  which  simply  demonstrates  the 


statistical  equilibration  obtained  by  multiple  measurement  in  the  presence  of  random 
error. 

The  1.3%  result  for  the  IDC  is  a  significant  improvement  over  the  3.3X  average 
absolute  error  obtained  using  the  same  raw  data  but  using  the  TDC.  The  most 
obvious  explanation  for  the  poorer  result  of  the  TDC  is  the  error  amplification 
suffered  as  a  consequence  of  the  higher  condition  number  of  ^Jj^with  the  TDC.  A 

rationalization,  untested  by  these  experiments,  is  that  the  TDC  uses  qc  to 
calculate  each  row  of  A^whereas  the  IDC  uses  each  measurement  to  calculate  no  more 

than  two  rows  of Since  q0  in  the  present  case  contains  values  nearest  the 
sensitivity  limit  of  the  instrument  (lowest  concentrations),  they  could  contain 
the  largest  errors  and  have  a  compounded  effect  on  the  final  results  with  the  TDC. 

The  overall  accuracy  of  the  PGSAM  lies  intermediate  between  the  TDC  and  the 
IDC.  Its  result  is  better  than  the  TDC,  most  probably  because  of  its  lower  error 
amplification  factor  [cond  £Al£>  -  1.0].  It  is  poorer  than  th'e  IDC,  most 

probably  because  of  experimental  considerations  explained  earlier. 

The  errors  obtained  by  the  IDC  using  the  second  group  of  wavelengths  are  sur¬ 
prisingly  low  considering  the  K  obtained  and  its  condition  number  (Table  V) . 

The  same  is  true  for  other,  even  less  selective,  wavelength  combinations  that 
allowed  concentrations  to  be  calculated  at  the  cost  of  analysis  accuracy. 

Finally,  as  a  direct  result  of  using  the  Househoulder  method,  a  more  detailed 
error  analysis  (see  Appendix)  was  obtained  and  is  given  in  Table  VII.  The  present¬ 
ation  of  these  results  (see  estimate  45,  Appendix)  allows  an  understanding  of  how 
each  source  of  error,  standard  addition  error  (£g)  and  measurement  error  (Ag), 
influences  the  determination  of  each  row  of  £  and  hence  the  final  concentrations. 

In  each  case,  the  generalized  condition  of^g  dominates  the  coefficient  of  the  err 
<£g)  term  to  give  that  term  a  greater  Influence  than  the  influence  of  the  err  (A^) 
term.  This  is  especially  true  of  the  TDC  method.  According  to  45  (Appendix), 

Table  VII  indicates  that  the  PGSAM  has  the  predicted  advantage  over  the  normal  GSAM 
even  using  the  IDC.  However,  for  reasons  given  earlier,  use  of  the  PGSAM  may  Indeed 


16 


lead  to  higher  values  for  err  (Aq)  and  perhaps  even  err  (Atj)  as  was  most  certainly 
found  to  be  true  here. 

CONCLUSION 

In  this  paper  the  theoretical  basis  for  the  constraints  on  multianalyte 
analysis  in  the  presence  of  sample  related  interferences  has  been  investigated. 

The  detailed  error  estimation  afforded  by  the  Househoulder  method  and  exploita¬ 
tion  of  the  properties  of  matrix  and  vector  norms  have  lead  to  improvements  to 
the  GSAM  as  published  in  an  earlier  paper  (2).  These  include: 

1)  A  new  definition  of  the  problem  to  be  solved  (the  IDC-GSAM)  and  a  scaling 
method  (equation  27)  together  yielding  a  vastly  reduced  error  propagation  for  the 
determination  of 

2)  A  new  experimental  design  (the  PGSAM)  with  desirable  error  propagation 
properties. 

3)  A  single  quantity,  the  cond  (K) ,  that  can  be  used  to,  quantitatively, 
compare  two  analytical  methods  designed  to  perform  the  same  multicomponent 
analysis. 

Finally,  a  more  quantitative  understanding  of  the  limitations  faced  by  the 

analytical  chemist  has  been  achieved.  Armed  with  a  collection  of  "fully  selective" 

sensors,  the  analyst  can  use  the  simple  standard  addition  method  to  perform  an 

analysis  for  several  analytes  on  a  complex  sample  with  absolutely  no  amplification 

of  the  measurement  uncertainty  as  it  is  propagated  through  the  mathematical  methods 

to  the  initial  analyte  concentrations.  Using  the  appropriate  variation  of  the 

Generalized  Standard  Addition  Method  the  analyst  can  now  perform  the  same  analysis 

with  "non-selective  sensors"  and  also  minimize  the  effect  of  error  amplification 

in  the  determination  of  the  linear  response  constants,  K,  and  the  initial  analyte 
•mounts,  n„,  or  concentrations,  cs.  Additionally,  with  an  approximation  (from 

analysis  of  a  previous  sample  or  standard)  of  the  selectivity  of  the  analytical 
system  as  expressed  by  K,  the  analyst  can  determine,  in  advance,  the  cost  of 
non-selectivity  in  units  of  error  amplification  limits.  In  some  cases  the  cost 
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might  be  high  enough  to  force  the  analyst  to  conclude  that  an  analysis  with  a 
desired  accuracy  may  be  futile. 
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APPENDIX 


Improving  Numerical  Stability 

Mathematically,  equation  10  represents  a  system  of  n  linear  equations  with  r 
right-hand  sides  Aq^,...,Aqr  (being  the  columns  of  AQ)  corresponding  to  the  r  un¬ 
known  colums  kj,...tk  of  jC.  As  the  number  of  additions  normally  exceeds  the 
number  of  analytes  and  sensors,  this  system,  being  overdetermined,  will  generally 
have  no  exact  solutions  since  the  right-hand  sides,  Aq^,...,Aq  ,  involve  randomly 
distributed  measurement  errors.  Therefore,  generalized  solutions,  k^,  yielding  an 
optimal  fit  to  equation  10,  in  a  certain  sense  to  be  defined,  must  be  found. 

In  the  following,  the  column  subscript  j  of  k^  and  q^  is  suppressed  fixing 


an 


arbitrary  index  l<J<r  to  give, 
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AM  -  Ac^,  ANelRn,r,  kelRr,  A^e*".  28 

For  any  vector  jceIR n,  its  Euclidian  norm,  ||jc  ,  is  defined  by 

II  n  ,  r  2,1/2  ,  T  ,1/2 

IML  =  (  l  V  *  <£x)  •  29] 

1  i-1  1 

A  vector  Je  will  be  called  a  generalized  (or  least  squares)  solution,  in  numerical 
analysis  jargon,  of  equation  28  if  it  minimizes  the  defect  of  equation  28,  i.e.,  if 
II  ANk  -  A^j|2  ”  II  AN*  -  A^l^  30) 


Minimizing  the  defect  of  a  linear  system  with  respect  to  the  Euclidian  norm  is  well 
known  as  linear  regression  and  leads  to  the  so-called  normal  equations, 

<£NT£N)k  -  31^ 


which  are 


derived  by  partially  differentiating  the  function 


-  II AB*  -  ^ll2  -  32) 

*  2*  ANTANx.  -  2xT ANTAq  +  AqTAc^ 

and  by  setting  the  partial  derivatives  to  zero.  The  normal  equations  have  a  unique 

solution  Jt  if  and  only  if  has  full  rank  r  in  which  case  is  formally  given  by 

-  ^ANTAN)_1ANTAq.  33) 

T  -IT 

The  matrix  (£N  £N)  ^  is  called  the  generalized  (or  Moore-Penrose)  inverse  of 


The  reader  should  keep  in  mind  that  the  definition  of  a  generalized  solution, 
k,  is  based  on  the  special  choice  of  the  Euclidian  norm  which  is  reasonable  in  the 


case  of 


normally  distributed ,  random  perturbations  of  the  right-hand  side 


entries  Aq.  (6).  In  the  case  of 


non-normally  distributed  errors,  the 


Lj-nortn 


H*  Hi  =  l  l*J 

1-1  1 


is  an  appropriate  function  to  be  minimized  (7),  i.e.. 


ll^Jk  -  A^Hj  -  min  ||£g&  -  A^. 
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Equation  28  might  also  be  poorly  conditioned  meaning  heuristically  that  a  small 
relative  error  of  the  right-hand  side,  produces  a  large  relative  error  of  the 
solution  vector  In  that  case  there  may  appear  large  negative  components  of 

k,  which  can  be  prevented  by  minimizing  equation  32  under  the  regularizing  con¬ 
ditions  >  0  (non-negative  least  squares  (ref.  8,  Chap.  23).  Also  a  combination 
of  the  L^-norm  and  the  non-negativity  conditions  might  be  reasonable. 

In  the  following,  consideration  is  restricted  to  the  simpler  problem  (equa¬ 
tion  30)  which  can  be  solved  by  a  very  effective  numerical  algorithm  allowing  a 
transparent  and  useful  error  analysis.  Therefore,  the  subscript  of  the  norm  is 
suppressed  writing  ||jc||  instead  of  HjjJL. 

At  first  sight,  the  most  obvious  way  to  compute  the  generalized  solution  of 

equation  30  is  to  solve  the  normal  equation  31.  In  fact,  if has  full  rank, 

T 

then  is  positive  definite,  and  almost  any  algorithm  appropriate  for  positive 

definite  linear  systems  (e.g.,  Cholesky's  algorithm)  may  be  applied.  From  the 
point  of  view  of  numerical  stability,  however,  it  is  well  known  (4,8)  that 
Househoulder' s  orthogonalization  is  highly  preferable  in  most  cases.  In  addition, 
Househoulder ' s  procedure  provides,  as  a  by-product,  some  useful  information  for 
an  a  posteriori  error  estimate  to  be  discussed  in  the  next  section.  A  short  out¬ 
line  of  the  algorithm  follovs: 

A  square  matrix  ^e!Rr,r  is  called  orthogonal ,  if  its  columns  span  an  ortho- 
IT  T 

normal  basis  of  1R  (i.e.,  11  *  i) •  The  product  of  orthogonal  matrices  is 
orthogonal.  Multiplication  of  a  vector  xeTRr  by  an  orthogonal  matrix  £  does  not 
change  its  Euclidian  length.  Therefore,  solving  the  minimization  problem  (equa¬ 
tion  30)  is  equivalent  to  solving 

II  ~  |^ll  •  «in  II  j*£Kx  -  PA^H 

xc»r 

for  any  orthogonal  matrix  P.  Mow,  a  matrix  P  of  the  form 


36) 
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T  T  r 

^-<I-2ww,w^-l,  we®  ,  37) 

y 

is  easily  shown  to  be  orthogonal  and  symmetric.  To  a  given  vector  ^e®  ,  ^  can 
be  chosen  In  a  way  that 

T  r 

•  Xe^,Xe®  ,  e^  ■  (1,0, . . .  ,0)  e®  .  38) 

Taking  the  first  column  of  for  in  the  first  step  and  writing  for  J>  the 
product  P^Jj  has  the  form 


P  Ajl 


39) 


the  stars  denoting  arbitrary  entries .  Applying  the  same  procedure  as  indicated 
above  (leaving  the  first  column  unaffected)  to  the  submatrix,  produces,  after 
r-1  multiplications  by  orthogonal  matrices  P^,...,Pr_1,  a  matrix  of  the  form 


where  ^  is  upper  triangular.  Now  equation  36  can  be  written  as 


Partitioning  the  vector  PAq  analogously  to  PAN  one  obtains  for  equation  41 


the  solution  of  which  is 

Jfe  -  ^  ■  ^1  43) 

As  is  triangular,  ^  can  be  calculated  explicitly  from  equation  43  provided^ 

has  full  rank  (in  which  case  £  is  nonsingular).  The  minimal  norm  of  the  defect 

according  to  equation  42  is 
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IlM  Ji  -  “  11^  II-  44) 

Since  equation  28  and  respectively  equation  30  must  be  solved  for  r  right-hand 
sides  Aq  , ...,Aa,  the  orthogonal  transformations  J)Aq^, .. jPAq  should  of  course  be 
performed  simultaneously  with  the  transformation  of  AN. 

Error  Estimation 

An  informative  error  bound  estimate  is  given  by  Reference  5,  p.  177.  It 
considers  matrix  errors  and  is  deduced  by  Taylor  expansion  neglecting  higher  order 
terms: 


2  II  M 

err^k)  <_  [cond^W)  +  condQj,)  ]  err  (AN)  + 


cond (W) 


ll  Aqll 

II  an|!  II ill 


II  an||  Ilk  II 

err  (Aq), 


where  the  relative  error  of  a  quantity,^,  err^) 


45) 


the  bar  indicating  the  perturbed  values.  £  and  h?  are  the  terms  appearing  in  the 
Househoulder  algorithm  (see  equations  43,  44).  Since  multiplication  by  an  ortho¬ 
gonal  matrix  leaves  the  norm  of  a  vector  invariant,  the  norm  of  AN  is  easily 
calculated  in  terms  of  W: 

ll Mil  =  maxIIMsll =  ®axIli,^ll e  max!lj&!l  Hljtll 

IIjeII  ml  llill-i  llill-i  46) 

As  W  is  a  small  triangular  matrix  the  computation  of  its  condition  and  its  norm  by 
determining  its  singular  values  is  easier.  Very  popular  (but  not  always  realistic) 
estimates  of  ||W  ||  and  cond  ^jj)  are  given  by  the  following  formulas: 

max  l w  | 

||W||  >  max  lwijl;  cond  ©  -  1,J  '  47) 

i,j  -l,...,r  min  |w  | 

i  11 


Thus  using  the  Househoulder  algorithm,  the  terms  appearing  on  the*  right  hand  side 
of  45  are  obtained  as  an  easily  calculable  byproduct  allowing  the  determin¬ 

ation  of  an  upper  bound  for  the  error  entering  the  columns  k^,...,k  of  the  £ 


matrix. 


TABLE  I 


Experimental  Design  and 

Responses  (AQ)  for  the  TDC-GSAM 


Ni 

Cu 

Co 

Cr2°7 

Aq(XA) 

Aq(X2) 

Aq(A3) 

Aq(X4) 

Aq(A5) 

1.00 

0.00 

0.00 

0.00 

0.0422 

0.491 

-0.010 

0.174 

0.058 

2.00 

0.00 

0.00 

0.00 

0.0822 

0.979 

-0.012 

0.345 

0.093 

3.00 

0.00 

0.00 

0.00 

0.0151 

1.48 

0.004 

0.532 

0.198 

4.00 

0.00 

0.00 

0.00 

0.209 

1.99 

0.018 

0.723 

0.266 

4.00 

2.00 

0.00 

0.00 

0.265 

1.98 

0.005 

1.11 

1.81 

4.00 

4.00 

0.00 

0.00 

0.326 

2.00 

0.002 

1.49 

3.28 

4.00 

6.00 

0.00 

0.00 

0.413 

2.00 

0.017 

1.89 

4.88 

4.00 

8.00 

0.00 

0.00 

0.507 

2.04 

0.044 

2.31 

6.44 

4.00 

8.00 

4.00 

0.00 

0.565 

2.16 

1.99 

2.41 

6.52 

4.00 

8.00 

8.00 

0.00 

0.502 

2.17 

3.85 

2.45 

6.39 

4.00 

8.00 

12.00 

0.00 

0.606 

2.32 

5.77 

2.59 

6.61 

4.00 

8.00 

16.00 

0.00 

0.733 

2.54 

7.80 

2.75 

6.70 

4.00 

8.00 

16.00 

8.00 

6.71 

4.22 

7.87 

2.71 

6.57 

4.00 

8.00 

16.00 

16.00 

12.6 

5.85 

7.90 

2.61 

6.18 

4.00 

8.00 

16.00 

24.00 

18.9 

7.70 

8.06 

2.65 

6.38 

4.00 

8.00 

16.00 

32.00 

25.4 

9.99 

8.45 

2.92 

6.58 
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TABLE  II 

Experimental  Design  ^N)  and 
Responses  (Ag)  for  the  IDC-GSAM 


Ni 

Cu 

Co 

Cr2°7 

Aq(Xx) 

Aq(X2) 

Aq(X3) 

Aq(X4) 

Aq(X5) 

1.00 

0.00 

0.00 

0.00 

0.042 

0.491 

-0.010 

0.174 

0.059 

1.00 

0.00 

0.00 

0.00 

0.040 

0.488 

-0.001 

0.172 

0.035 

1.00 

0.00 

0.00 

0.00 

0.069 

0.497 

0.016 

0.186 

0.105 

1.00 

0.00 

0.00 

0.00 

0.058 

0.516 

0.013 

0.191 

0.068 

0.00 

2.00 

0.00 

0.00 

0.056 

-0.007 

-0.012 

0.382 

1.54 

0.00 

2.00 

0.00 

0.00 

0.060 

0.018 

-0.003 

0.388 

1.47 

0.00 

2.00 

0.00 

0.00 

0.087 

-0.008 

0.015 

0.399 

1.60 

0.00 

2.00 

0.00 

0.00 

0.094 

0.048 

0.027 

0.419 

1.56 

0.00 

0.00 

4.00 

0.00 

0.058 

0.115 

1.95 

0.103 

0.079 

0.00 

0.00 

4.00 

0.00 

-0.063 

0.015 

1.86 

0.034 

-0.129 

0.00 

0.00 

4.00 

0.00 

0.104 

0.152 

1.92 

0.144 

0.223 

0.00 

0.00 

4.00 

0.00 

0.128 

0.218 

2.03 

0.154 

0.089 

0.00 

0.00 

0.00 

8.00 

5.98 

1.67 

0.075 

-0.039 

-0.125 

0.00 

0.00 

0.00 

8.00 

5.93 

1.63 

0.026 

-0.098 

-0.389 

0.00 

0.00 

0.00 

8.00 

6.27 

1.85 

0.165 

0.036 

0.192 

0.00 

0.00 

0.00 

8.00 

6.50 

2.30 

0.389 

0.275 

0.202 

TABLE  III 


Experimental  Design  and 

Responses  (AQ)  for  the  PGSAM  (TDC) 


Ni 

Cu 

Co 

Cr2°7 

Aq(X1) 

Aq(X2) 

Aq(X3) 

Aq(X4) 

Aq(X5) 

1.00 

0.00 

0.00 

0.00 

0.078 

0.513 

0.026 

0.193 

0.103 

2.00 

0.00 

0.00 

0.00 

0.100 

1.00 

0.010 

0.363 

0.129 

3.00 

0.00 

0.00 

0.00 

0.172 

1.52 

0.029 

0.552 

0.292 

4.00 

0.00 

0.00 

0.00 

0.220 

2.02 

0.025 

0.732 

0.286 

0.00 

1.00 

0.00 

0.00 

0.031 

-0.006 

0.001 

0.200 

0.799 

0.00 

2.00 

0.00 

0.00 

0.072 

0.013 

0.010 

0.393 

1.52 

0.00 

3.00 

0.00 

0.00 

0.108 

0.021 

0.019 

0.596 

2.31 

0.00 

4.00 

0.00 

0.00 

0.166 

0.039 

0.028 

0.803 

3.09 

0.00 

0.00 

1.00 

0.00 

0.010 

0.023 

0.495 

0.028 

0.035 

0.00 

0.00 

2.00 

0.00 

0.036 

0.050 

0.962 

0.047 

0.019 

0.00 

0.00 

3.00 

0.00 

0.020 

0.080 

1.45 

0.078 

0.058 

0.00 

0.00 

4.00 

0.00 

0.019 

0.104 

1.94 

0.100 

0.066 

0.00 

0.00 

0.00 

1.00 

0.718 

0.197 

0.011 

-0.002 

0.005 

0.00 

0.00 

0.00 

2.00 

1.51 

0.435 

0.026 

-0.004 

-0.030 

0.00 

0.00 

0.00 

3.00 

2.22 

0.642 

0.035 

-0.003 

-0.009 

0.00 

0.00 

0.00 

4.00 

3.07 

0.891 

0.048 

-0.005 

-0.018 

TABLE  IV 

K  Calculated  by  IDC-GSAM  Using  Wavelengths 


V 

A2,  *3»  ^5* 

Cond  ££) 

-  1.93 

^2 

h. 

Cr 

0.76 

0.23 

0.02 

-0.01 

Ni 

0.05 

0.49 

0.00 

0.06 

Co 

0.08 

0.02 

0.48 

0.01 

Cu 

0.04 

0.01 

0.01 

0.77 
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TABLE  V 

K  Calculated  by  the  IDC-GSAM  Using  Wavelengths 


^1*  X^» 

X^,  Xj. 

Cond  QC) 

-  5.24 

h 

^3 

Cr 

0.77 

0.02 

0.01 

-0.01 

Co 

0.01 

0.48 

0.03 

0.02 

Ni 

0.05 

0.01 

0.18 

0.07 

Cu 

0.04 

0.01 

0.20 

0.77 

Generalized  Cond  ^AN) 


T 

Eigenvalues  of  K  K  1 

2 

3 

4 


Cond  (K) 


Calculated  Initial  Units 


Ni 

(True  -1.0) 

Cu 

(True  -  1.0) 

Co 

(True  «  1.0) 

Cr 

(True  -  0.0) 

TABLE  VI 
Final  Results 


^2*  ^3*  ^5 

TDC-GSAM  IDC-GSAM 


3.42 


0.701 

0.586 

0.231 

0.195 


1.90 


1.034 

1.015 

1.032 

0.018 


Average  Absolute  X  Error  *  3.30 


Standard  Error  Estimate  **  I  3.0 


1.0 


0.723 

0.585 

0.235 

0.194 


1.93 


1.018 

1.002 

1.014 

0.006 


1.33 


PGSAM 


1.0 


0.700 

0.579 

0.236 

0.198 


1.88 


0.995 

0.962 

0.989 

0.015 


2.30 


2.5 


*1»  *4» 

IDC-GSAV 


1.0 


0.666 

0.582 

0.234 

0.024 


5.24 


1.032 

1.001 

1.014 

0.005 


1.73 


2. 


_  _  .  „„  „  £  |  True-Calculated  | 

*  Average  Absolute  X  Error  -  100.0  x  £(lnltlai  amounts) 

.  .  _  .  _  _  .  Fz  (True-Calculated) 2  I  *^2 

**  Standard  Error  Estimate  -  |Z(lnitial  amounTs)2 


x  100.0 
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TABLE  VII 

Error  Analysis  According  Co  Equation  45 

Method  Wavelengths  Coefficients  of  err  (AN)*  Coefficients  of  err  (Aq)** 


TDC-GSAM 

v 

V 

X3* 

X5 

3.42,  3.42, 

3.42,  3.42 

0.54, 

0.90, 

0.70, 

0.78 

IDC-GSAM 

V 

X2’ 

X3* 

X5 

1.02,  1.04, 

1.03,  1.04 

0.51, 

0.51, 

0.50, 

0.51 

PGSAM 

xl* 

X2* 

X3* 

X5 

1.0,  1.0,  1 

.0,  1.0 

0.18, 

0.18, 

0.18, 

0.18 

IDC-GSAM 

» 

1.02,  1.03, 

1.04,  1.04 

0.50, 

0.50, 

0.50, 

0.50 

2  II  h2  11 

*  [cond  (W)  +  cond  (W)  - - -  ] 

II  AN||  II  k|| 

II  Aq  H 

||AN||  ||  k|| 


**  [cond  (W) 


ta* 


. 
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